Skip to content

plt.pcolormesh(x_mesh/1e3-x0, y_mesh/1e3, marmousi_inv)

plt.figure(figsize=(12,8)) plt.pcolormesh(x_mesh/1e3-x0, y_mesh/1e3, marmousi_std/marmousi_inv * 100, vmin=0, vmax=np.max(marmousi_std/marmousi_inv*100)/1.3, rasterized=True) # plt.pcolormesh(x_mesh/1e3-x0, y_mesh/1e3, marmousi_inv) # plt.plot(x[receiver[0][0]]/1e3-x0, y[receiver[0][1]]/1e3, 'wv', alpha=0.5) # for i in range(len(source)): # plt.plot(x[source[i][0][0][0]]/1e3-x0, y[source[i][0][0][1]]/1e3, 'r*') plt.xlabel("X (km)") plt.ylabel("Z (km)") plt.gca().invert_yaxis() plt.axis('scaled') cb = plt.colorbar(shrink=0.5) # cb.ax.set_title('std($V_p$) (km/s)', fontsize=14) cb.ax.set_title('%', fontsize=14) plt.savefig("marmoursi-UQ-std0-norm.png", bbox_inches='tight', dpi=600) plt.savefig("marmoursi-UQ-std0-norm.pdf", bbox_inches='tight', dpi=600) plt.show()